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Disappearance of a stacking fault in the hard-sphere crystal under gravity, such as reported by Zhu et al. 
[Nature 387 (1997) 883], has successfully been demonstrated by Monte Carlo simulations. We previously found 
that a less ordered (or defective) crystal formed above a bottom ordered crystal under stepwise controlled 
gravity [Mori et al. J. Chem. Phys. 124 (2006) 174507]. A defect in the upper defective region has been 
identified with a stacking fault for the (001) growth. We have looked at the shrinking of a stacking fault 
mediated by the motion of the Shockley partial dislocation; the Shockley partial dislocation terminating the 
lower end of the stacking fault glides. In addition, the presence of crystal strain, which cooperates with gravity 
to reduce stacking faults, has been observed. 
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1 Introduction 

The existence of a crystal-fluid phase transition 
in the hard-sphere (HS) system was observed in 
1957 using both molecular dynamics (MD) [T] and 
Monte Carlo (MC) simulation techniques [2]. In 
the HS system the phase behaviour is governed by 
the particle number density only; the disordered 
fluid phase crystallizes into a face-centred cubic 
(fee) crystal over a coexistence interval 0.494 < (f) 
< 0.545 fSH, where (t)={Tr/6)a^{N/V) is the volume 
fraction of the HSs. Here, a is the HS diameter, N 
the number of particles, and V the volume of the 
system. The HS system is often utilized as a mod- 
els for colloids [5]. In order for colloidal crystals 
to have application as photonic crystals O [3 [8], 
it is important that the number of defects, such as 
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II By adjusting the lateral box size so as to remove the 
stress of the crystal the coexistence region has been corrected 
slightly, i.e., 0.491 < <j) < 0.542 [1]. However, the result has 
not been affected; the stress cannot be removed entirely be- 
cause the system is nonuniform and, in addition, the crystal 
can rotate about z axis 



stacking faults, be reduced. However, the number 
of stacking faults present in a HS crystal is diffi- 
cult to reduce, because the stacking disorder does 
not change the density. Indeed, the entropy differ- 
ence between fee and hexagonal close-packed (hep) 
HS crystals has been calculated [TOl [TTl [II E] 
and is, at most, of the order of 10~'^A;b per par- 
ticle, where /cb is Boltzmann's constant. On the 
other hand, the density gradients induced by grav- 
ity have been observed to enhance colloidal crys- 
tallization. For example, HS-like suspensions, such 
as poly(methylmethacrylate) (PMMA) microparti- 
cles dispersed in a hydorcarbon medium as well as 
charge stabilized colloidal suspensions, have been 
seen to crystallize by sedimentation 114j. 

Interestingly, the 1997 experiment of 
Zhu et al. [15] suggests that stacking faults in 
HS suspensions were controlled by gravity. In 
this paper, we focus on the mechanism of the 
reduction of the stacking faults in colloidal crystals 
due to gravity. The conclusion of Zhu et al. is 
as follows: Whereas the resulting HS crystalline 
structure under microgravity on the Space Shuttle 
was observed to be random hexagonal-close-packed 
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(rhcp) [15] ■ on Earth a mixture fee and rhep ean 
be obtained [16] . The resuhs of the experiments of 
Kegel and Dhont [l7] support this trend, although 
the final structure was the faulted- twinned fee [TS]. 
The metastability of the rhcp structure and sta- 
bility of the fee structure have been supported by 
numbers of experimental results [191 IlQl [111 [22] . In 
addition to the reduction of the stacking disorder, 
these results also indicate that some stacking 
disorders remain at a late stage. The difference 
between the states at a late stage implies that the 
states observed are not in true equilibrium. Thus, 
the disappearance of defects in most cases must be 
a process between two metastable states. Though 
what this paper deals may be mere one of such 
processes, it should provide an insight into con- 
trolling defects. A recent remarkable experimental 
result, which supports the variety of the metastable 
states, is an small-angle X-ray diffraction study of 
silica colloids by Dolbnya et al. ^23j. Their samples 
were aged for 1.5-2 years, which is in contract to 
those in previous studies: the times involved in 
the experiments of Zhu et al. were limited by the 
flight time of the Space Shuttle, although details 
were not described, and the aging times in the 
experiments of Kegel and Dhont were, at most, a 
few tens of days. It should be pointed out that [20] 
have presented a result, not reported in [15], that 
even under microgravity the fee structure develops 
after thousands of seconds. Dolbnya et al. [23] 
have found a scattering pattern indicating a mostly 
complete fee. Similar results have been obtained 
for a PMMA dispersion by a light scattering by 
Martelozzo et al. [H]. In opposition to the con- 
clusion of Dolbnya et al. that melting of the rhcp 
structure followed by re-growth of fee crystal is 
likely, we stress that this scenario is only one of the 
many processes possible. Despite that pre-melting 
near many types of defects, including the Shockley 
partial dislocation, has been observed [24^ and 
supports this scenario, our opinion is that gradual 
disappearance of a defect, mediated by the motion 
of related defects, such as the glide of dislocation, 
cannot be ruled out. 

The growth of a {lll}-oriented crystal occurs 
more often on a flat, bottom wall. This phe- 
nomenon can be understood through the orien- 
tational dependence of the interfacial free energy, 
which for an fee HS crystal against a hard flat 
wall is lowest for a {111} crystal orientation at 
the wall, compared with the other two common 
low-index orientations: {100}, and {110}) [25]. In 
the {111} growth, however, stacking disorder can- 



not be avoided for the reason discussed above; 
that is, the density is the same for any stack- 
ing sequences. For {100} growth, on the other 
hand, the stacking sequence is unique. In a tech- 
nique sometimes referred to as colloidal epitaxy, 
van Blaaderen et al. [26], used a template to force 
the crystal growth in the <001> direction. A pat- 
terned substrate [221 [571 UHl US] or a substrate 
with a structure possessing the symmetry same as 
the {100} face [SO] |3T1 [32] must provide a stress 
as a driving force for the {100} growth. Despite 
this bias toward {100} growth, Schall et al. [33] 
used confocal microscopy to observe a stacking fault 
running along {111} in the colloidal epitaxy of a 
PMMA suspension. This stacking fault was found 
to be terminated by the Sockley partial dislocation. 
Schall et al. interpret this as a misfit dislocation; 
thus, this structure is in equilibrium in the sense 
that the strain energy is minimized. In this pa- 
per, we stress that the same structure, which is in 
a metastable state, does disappear. One can say 
that the {001} growth has an additional advantage 
that even if some excess stacking faults exist in a 
stage, some of them may disappear in a later stage 
through the mechanism we shall present in this pa- 
per. Again, the purpose of the paper is to provide 
a mechanism for the gradual disappearance of a de- 
fect. A detailed statistical analysis of this mecha- 
nism will be left to future research. 

2 Simulation 

In previous MC simulations [34], we observed the 
formation of a highly defective crystalline region 
above a well-ordered crystal at the bottom of a HS 
system under gravity. The structure of the highly 
defective crystalline region depends on the growth 
direction: For {100} growth the highly defective 
region was observed to include several kinds of de- 
fects, whereas for the {111} growth the defective 
region was smectic-like in appearance, and, upon 
closer examination, was observed to consist of layers 
with significant crystallinity. For the {100} growth, 
we used a system size of = 1664 particles con- 
tained in a box of dimensions Lx=Ly^6.27a and 
^^==49.2(7. For the {111} growth, N = 3744 parti- 
cles were used with Lx=Ly=9.4:0a and Lz—GO.Ocr. 
Periodic boundary conditions were used in the x 
and y directions and well-separated hard walls were 
placed at at z — (top) and z — (bottom). 
We speculate that the observed {001} growth re- 
sults from stress due to the small system size and 
square symmetry in the horizontal (x and y) direc- 
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tions. This stress must play an similar role in the 
colloidal epitaxy [26l[2Il[28l[29l[30l[3ll[32l[33]. At 

least, the differences in the diagonal components of 
the strain from one another is of the same symme- 
try because, whereas the horizontal lattice spacing 
is fixed, the vertical one can vary |35j . 

In this study [34] the magnitude of the grav- 
itational force was increased stepwise; the gravi- 
tational constant g*='mga/k^T was increased by 
an increment A(7*=0.1 every 2x10^ MC cycle 
(MCC) (At=2xlO^MCC) up to a maximum value 
of g* — l.b. Here, m is the mass of a particle, T the 
temperature, and one MCC is defined by a set of N 
position moves. We remark that through this grad- 
ual stepwise increase of g* we can successfully avoid 
metastable polycrystalline (we have suffered from 
these phenomena in an early stage of the study [36]) 
and glassy states. The effect that the gravity slows 
down the crystallization at high density has been 
observed by a MD simulation I37j . in agreement 
with [36,. Also we note that the stepwise control 
can be efficiently realized by use of a centrifugation 
rotator [38] [3^. In earlier work, we reported that 
the formation of a less defective crystalline region 
started around 5'*=0.9; up to g*=0.8 the bottom of 
the system was occupied by a defective crystal [34] . 
We report here the behaviour of the system, which 
exhibited the (001) growth, focusing on the stacking 
faults. 

3 Shrinking stacking fault 

To illustrate the evolution of the system configura- 
tion for g* — 0.9, we show in Fig. [1] (a)-(e), three- 
dimensional snapshots at 1.90, 1.92, 1.94, 1.96, and 
1.98x lO^MCC, respectively. The regions where dis- 
order in stacking sequence occurs are marked with 
a red line except for Fig. [T] (a) ; we shall identify this 
disorder with an intrinsic stacking fault later. Based 
on these snapshots, we can conclude that shrinkage 
of an intrinsic stacking fault occurs in this system. 
In Fig. [1] (f) we plot the center of the gravity, zq, 
which moves to lower positions as the simulation 
proceeds (though the '4ime" average is not taken, 
the data is sufhcient to show the trend) . While the 
single stacking fault is shrinking, the velocity of de- 
creasing of zq is small. On the other hand, the 
sinking velocity is high when the conversion rate 
from the disordered to ordered states is high. In 
addition, we observe that the motion is not simply 
associated with the increasing height of the deficit 
region (as we shall see the lower end of an intrinsic 
stacking fault accompanies a particle deficiency as 



shown in Fig. [3]) in a viscous medium; the motion is 
better explained in terms of the elastic interaction 
between the defect and the strain field. In some 
cases buoyancy and elastic effects cooperate with 
each other to accelerate the motion. We have ob- 
served that the conversion of the defective crystal 
into an ordered crystal under gravity is caused by 
disappearance of stacking faults. In other words, we 
have successfully demonstrated the role of gravity 
reducing stacking faults. One point that should be 
noted is that the observations have been made over 
a relatively short period, along with a short simula- 
tion duration. We again note that there are many 
possible processes involved in the disappearance of 
defects and what we are observing here are those 
processes that take place at relatively short time 
scales. In addition, because the magnitude of At 
is not large, it follows that the stepwise control of 
the centrifugation rate is promising in a sense that 
the mechanism described in this paper must be of 
realizable time scale. 

4 Shockley partial dislocation 

Evolution of the particle configuration within the 
18th to 21st layers from the bottom of the system 
is shown in Fig. [2] The projections of the configura- 
tions onto xy plane correspond approximately to a 
square lattice, indicating that the lattice planes are 
fee (001) planes with the [110] and [110] directed 
along the x and y directions. The shrinkage of the 
defect in Fig. [1] corresponds to the disappearance 
of regions where the square lattice is shifted by a 
half of its lattice constant, a/v^, where a is the 
fee lattice constant (a/V2 is equal to the nearest- 
neighbour interparticle separation for an fee crys- 
tal). In particular, disappearance of the defect is 
completed in the 18th layer at 1.94 xlO^th MCC 
and in the 19th layer at 1.96 xlO^th MCC. The 
magnitude of this shift is equal to the length of 
projection of the Burgers vectors of the Shockley 
partial dislocations onto the direction of the Burg- 
ers vector of the corresponding perfect dislocation. 
The translation of a fee (111) plane along [101] by 
one lattice constant, a/\/2, of the triangular lat- 
tice of the (111) plane is decomposed into two mo- 
tions: one of which translates particles in A-type 
sites onto B-type sites and the other in B-type to 
adjacent A-type [40]. Here, A and B refer to two 
of the three possible layer positions in projection 
onto the (111) plane (C refers to the last one). The 
translation vector b, which is the Burgers vector 
of a perfect dislocation, for one interparticle spac- 
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Figure 1: Figure[TJ Evolution of particle configuration in three-dimensional view and the center of gravity. 
We show views from -x direction at (a) 1.90, (b) 1.92, (c) 1.94, (d) 1.96, and (e) 1.98 xlO^th MCC. The 
red lines indicate a stacking fault, which is shrinking as simulation proceeds [The line is not drawn in (a) 
because the structure is complicated.] (f) the center of gravity, zq, is plotted. 



ing along <110> direction is decomposed into b + 
6^^, where and 6^^ are the Burgers vectors of the 
Shockley partial dislocations. Translation by b or 
b'^ also gives rise to the shift of lattice planes in 
the vertical direction that is shown in Fig. [3] - the 
blue (dashed) lines are shifted upward with respect 
to the green (dotted) lines. In Fig. [1] we examine 
a surface of the simulation box from —x direction. 
The fee unit cell is, thus, as illustrated in Fig. IH 
The Burgers vectors are, respectively, b =(1/2)[110] 
(=ai/2+a2/2 with ai, 02, and 03 denoting the lat- 
tice vectors), b^ = (1/6)[211], and 6" = (1/6)[12T]; 
6^ which makes a shift of magnitude a/2\/2 in the 
X direction in Fig. [2] and a shift of magnitude a/6 
in the vertical direction in Fig. [31 is drawn by an 
arrow in Fig. [H The difference in the magnitudes 
of the shifts in horizontal and vertical directions in 
Fig. [3] is affected by the differences in lattice con- 
stants in the vertical and horizontal directions due 
to the compression of the fee unit cell in the ver- 



tical direction |35j . In this way, we show that the 
region marked by blue (dashed) lines in Fig. [3] slides 
along [211] direction. This operation yields an in- 
trinsic stacking fault, i.e., stacking sequence such 
as ABCBCA- • • is made from ABCABC- • •. In con- 
clusion, an intrinsic stacking fault shrank through 
a glide of the Shockley partial dislocation. We note 
that the Frank partial dislocations, whose Burgers 
vector in the present case was (1/3) [111], were not 
observed; if this had occurred the number of planes 
along [111] would have varied across the Frank par- 
tial dislocation. 

In addition to the stacking fault, we find deforma- 
tion of the crystal; the crystal lattice in the defec- 
tive region was twisted with respect to the bottom- 
ordered lattice as shown in Fig. [2] and the lattice 
planes were inclined with respect to the bottom- 
ordered crystal as shown in Fig. [31 Twist and tilt 
boundaries are anticipated along the horizontal di- 
rection. Comparing the 19th and 20th layers in 
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Figure 2: Figure [2j Evolution of particle configura- 
tion of several crystalline layers. We plot particle 
positions in 18tli to 21st layers from the bottom of 
the system at (a) 1.92, (b) 1.94, and (c) 1.96 xlO^ 
th MCC. The z coodinate of each layer can be taken 
from Fig. [31 



Fig- [21(c) with those in Figs. [21 (a) and (b) the twist 
is found to disappear as the stacking fault disap- 
pears. Also, we see the simultaneous disappearance 
of the inclination. Taking a close look at Fig. [31 
we see that, due to the inclination of the lattice 
planes, the right side and left sides of a lattice plane 
appear to match up with each other across the pe- 
riodic boundary. This means that, for instance, the 
B-typc (111) plane gradually changes to the A- type 
(111) plane. It follows that strain energy is released 
during the process of stacking fault disappearance. 
Of course, there is a density increase correspond- 
ing to the translation of lattice planes downward. 
It is conjectured that the cooperation of those ef- 
fects makes the disappearance of stacking faults at a 
condition with g*—0.9] this condition can be rewrit- 



T 



t # 

• • t 





y/CT 



Figure 3: Figure [3l Side view of crystalline layers: 
18th to 22nd layers are from the bottom. The re- 
gion -2.92 < x/a < -0.92 at 1.94 xlO^th MCC 
is viewed from —x direction; a few layers from the 
surface are shown to avoid complication. Dashed 
(blue) and dotted (green) lines, drawn for guide for 
viewing, indicate the slide of dashed (blue) region 
occurs along the thick (red) line. 



ten so that the gravitational length l^—k^T /nig is 
slightly greater than the HS diameter tr. 

We should add a remark with regard to the strain. 
The interlayer spacing varies along the vertical di- 
rection whereas the interparticle separation in a 
layer does not vary appreciably [35j , so non-uniform 
uniaxal compression is observed. However, the re- 
lief of twist and tilt must lower the system free 
energy irrespective of the existence of non-uniform 
uniaxal compression. In our opinion, the defects 
in colloidal crystals can be controlled through the 
mechanisms found by this MC simulation. We note 
the effect of system size on the inclination of a 
layer; in contrast to the twists, stress due to the in- 
clination may decrease with increasing system size 
as the inclination angle becomes small or matches 
that corresponding to the periodicity of the "dislo- 
cation" array that minimizes the strain energy. 



5 Conclusion and remarks 

We have successfully demonstrated the disappear- 
ance of a stacking fault for (001) growth in a hard- 
sphere fee crystal under gravity In other words, we 
have given an answer to the speculation that stack- 
ing faults are reduced under gravity |15j . While an 
intrinsic stacking fault is terminated by the Shock- 
ley partial dislocation at an end, traversing the sys- 
tem horizontally the translation gradually disap- 
pears in a plane through the gradual variation of 
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—mgp{z) with p{z) being the particle number den- 
sity)]. We present the result in this style, with 
the strain energy calculation being left at the fu- 
ture research, is because the strain illustrated in the 
present figures evidently contribute dominantly; the 
energy of such strain field logarithmically depends 
on the linear dimension over which the field extends. 
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Figure 4: Figure 2] Schematic figure of translation 
corresponding to the Shockley partial dislocation. 
The (111) plane is marked (painted in red). The 
Burgers vector of the Shockley partial dislocation, 
b^= (1/6) [211] = ai/3 +02/6 +03/6, is shown by 
an arrow. Here, lattice vectors are denoted as ai, 
a2, and a^. We see that 6^ makes shifts of mag- 
nitude a/2-\/2, a/6\/2, and a/6, respectively, along 
[110], [110], and [001]. 

the magnitude of inclination of a (001) layer. In 
addition, we observe the existence of twists of por- 
tions of a (001) plane with respect to the bottom- 
ordered crystal. Those crystal deformations that 
disappear during shrinkage of a stacking fault do 
so via a glide of the Shockley partial dislocation. 
The disappearances of twist and tilt boundaries as 
well as of the deformations are considered to lower 
the system free energy. The density increase, which 
accompanies the decrease of the gravitational en- 
ergy during the shrinkage of stacking faults, has a 
tendency to cooperate with those phenomena. 

In this paper, we have provided some intuition 
into the gravity-induced disappearance of stacking 
faults in a hard-sphere crystal, while avoiding a de- 
tailed statistical analyses employing strain-energy 
fomulae and numerics. A strain-energy formula- 
tion for the present defect configuration, such as 
described in [40], is in progress [41]. On the other 
hand, some of parameters necessary for accurate 
numerical evaluation are lacking [namely, the dis- 
location core energy and the density dependence 
of the stacking fault free energy in the strongly 
inhonogenous case, such as in the high-g* regime 
(note that the pressure gradient, dP/dz, equals 

** The inerfacial free energy of stacking fault calculated 
near the crystal-fuid phase transition 1121 is of order 10"'' 
fcgT per unit area, so the interfacial energy contributes less 
as compared with the strain energy. 



References 

[1] B. J. Alder and T. E. Wainwright, 
J. Chem. Phys. 27, 1208 (1957). 

[2] W. W. Wood and J. D. Jacobson, 
J. Chem. Phys. 27, 1207 (1957). 

[3] W. G. Hoover and F. H. Ree, J. Chem. Phys. 
49, 3609 (1968). 

[4] R. L. Davidchack and B. B. Laird, 
J. Chem. Phys. 108, 9452 (1998). 

[5] P. N. Pusey, in Liquids, Freezing and 
Glass Transition, edited by J. P. Henderson, 
D. Levesque, J. Zinn- Justin, (Elsevier, Ams- 
terdam, 1991). 

[6] K. Ohtaka, Phys. Rev. B 19, 5057 (1979). 

[7] E. Yablonovitch, Phys. Rev. Lett. 58, 2059 
(1987). 

[8] S. John, Phys. Rev. Lett. 58, 2486 (1987). 

[9] D. Frenkel and A. J. C. Ladd, J. Chem. Phys. 
81, 3188 (1984). 

[10] L. V. Woodcock, Nature 385, 131 (1997); 
P. G. Bolhuis, D. Frenkel, S.-C. Mau, and 
D. A. Huse. ibid 338, 235 (1997) . 

[11] A. D. Bruce, N. B. Wilding, G. J. Ackerson, 
Phys. Rev. Lett. 79, 3002 (1997). 

[12] S. Pronk and D, Frenkel, J. Chem. Phys. 110, 
4589 (1999). 

[13] S.-C. Mau and D. A. Huse, Phys. Rev. E 59, 
4396 (1999). 

[14] K. Davis, W. B. Russel, and W. J. Glantsching, 
Science 245, 507 (1989). 



A. Mori et al. 



7 



[15] J. Zhu, M. Li, R. Rogers, W. Meyer, R. H. Ot- 
tewill, STS-73 Space Shuttle Crew, W. B. Rus- 
scl, and R. M. Chaikin, Nature 387, 883 
(1997). 

[16] R N. Pusey, W. van Megen, R. Bartlett, 
B. J. Ackerson, J. G. Rarity, and S. M. Un- 
derwood, Phys. Rev. Lett. 63, 2753 (1989). 

[17] W, K. Kegel and J. K. G. Dhont, 
J. Chem. Phys. 112, 3431 (2000). 

[18] C. Dux and H. Versmold, Phys. Rev. Lett. 78, 
1811 (1997). 

[19] M. S. Elliot, B. T. F. Bristol, and 
W. C. K. Poon, Physica A 235, 216 (1997). 

[20] Z. Cheng, J. Zhu, W. R. Russel, W. V. Meyer, 
and P. M. Chaikin, Appl. Opt. 40, 4146 (2001). 

[21] V. C. Martelozzo, A. B. Schofield, 
W. C. K. Poon, and P. N. Pusey, Phys. Rev. E 
66, 021408 (2002). 

[22] J. P. Hoogcnboom, D. Derks, P. Vergeer, and 
A. van Blaaderen, J. Chem. Phys. 117, 11320 
(2002). 

[23] I. P. Dolbnya, A. V. Petukov, D. G. A. L.Aarts, 
G. J. Vroege, and H. N. W. Lekkerkerker Eu- 
rophys. Lett. 72, 962 (2005). 

[24] A. M. Alsayed, M. F. Islam, J. Zhang, 
P. J. CoUings, and A. G. Yodh, Science 309, 

1207 (2005). 

[25] M. Hcni and H. Lowen, Phys. Rev. E 60, 7057 

(1999). 

[26] A. van Blaaderen, R. Ruel, and R. Wiltzius, 
Nature 385, 321 (1997). 

[27] K.-h Lin, J. C. Crocker, V. Prasad, 
A. Schofield, D. A. Weitz, T. C. Luben- 
sky, and A. G. Yodh, Phys. Rev. Lett. 85, 
1770 (2000). 

[28] D. K. Yi, E. M. Sco, and D.-Y. Kim, 
Appl. Phys. Lett. 80, 225 (2002). 

[29] J. P. Hoogcnboom, A. Ycthiral, 
A. K. van Langen-Suuring, J. Romijn, 
and A. van Blaaderen, Phys. Rev. Lett. 89, 
256104 (2002). 

[30] Y. Yin and Y. Xia, Avd. Mater. 14, 605 (2002). 



[31] J. Zhang, A. Alsayed, K. H. Lin, S. Sanyal, 
F. Zhang, W.-J. Pao, V. S. K. Balagu- 

rusamy, P. A. Hcincy, and A. G. Yodh, 
Appl. Phys. Lett. 81, 3176 (2002). 

[32] S. Matsuo, T. Fujine, K. Fukuda, 
S. Joudokazis, and H. Misawa, 
Appl. Phys. Lett. 82, 4283 (2003). 

[33] P. SchaU, I. Cohen, D. A. Weits, and 
F. Spaepen, Science 305, 1944 (2004). 

[34] A. Mori, S.-i. Yanagiya, Y. Suzuki, T. Sawada, 
and K. Ito, J. Chem. Phys. 124, 174507 (2006). 

[35] A. Mori, S.-i. Yanagiya, Y. Suzuki, T. Sawada, 
and K. Ito, Sci. Technol. Adv. Mater. 7, 296 
(2006). 

[36] S.-i. Yanagiya, A. Mori, Y. Suzuki, Y. Miyoshi, 
M. Kasuga, T. Sawada, K. Ito, and T. Inoue, 
Jpn. J. Appl. Phys., Part 1, 44, 5113 (2005). 

[37] I. Velkov, M. Cieplak, J. Koplik, and R. Ba- 
navar, Phys. Rev. E 66, 061401 (2002). 

[38] M. Megens, C. M. van Kats, P. Bosecke, and 
W. L. Vos, J. Appl. Cryst. 30, 641 (1997). 

[39] B. J. Ackerson, S. E. Paulin, B. John- 
son, W. van Megen, and S. Underwood, 
Phys. Rev. E 59, 6903 (1999). 

[40] J. H. Hirth and J. Lothe, Theory of Disloca- 
tion, 2nd. ed. (Krieger, Malaba, 1982). 

[41] A. Mori and Y. Suzuki, submitted. 



